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Production of electroweak bosons in e^e annihilation at high energies 
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Production of electroweak bosons in e^e~ annihilation into quarks and into leptons at energies 
much greater than 100 Gev is considered. We account for double-logarithmic contributions to all 
orders in electroweak couplings. It is assumed that the bosons are emitted in the multi-Regge 
kinematics. The explicit expressions for the scattering amplitudes of the process are obtained. It is 
shown that the cross sections of the photon and Z production have the identical energy dependence 
and asymptotically their ratio depends only on the Weinberg angle whereas the energy dependence 
of the cross section of the W production is suppressed by factor s~°''^ compared to them. 

PACS numbers: 12.38.Cy 



I. INTRODUCTION 



in the double- logarithmic approximation (DLA) was considered first in Ref . 0| . In this work 
y—i • it was shown that when the total energy of the annihilation is high enough, the most sizable radiative QED corrections 
J> ' to e+e~ annihilation into are double-logarithmic (DL). These corrections were calculated in Ref. Q] to all orders 

in a. The DL contributions to this process appear when the final /i"*"/^ -pair is produced in the Regge kinematics, i.e. 

■ when the muons move (in cmf) closely to the initial e+e^ -beam direction. According to the terminology introduced in 
1 ' , Ref. [0, the process where /x"*" ixioves in the (e~) -direction is called forward (backward) annihilation. Generalisation 

■ of these results to QCD (the forward and backward annihilation of quarks into quarks of other flavours) and to the EW 
, theory (the backward annihilation of the left handed leptons into the right handed leptons) was obtained in Ref. 

' and Ref. respectively. The forward and backward annihilation of e^e~ into quarks, all chiralities accounted for, was 
considered recently in Ref. 0|. One of the features obtained in Refs is that the forward scattering amplitudes 

in DLA are greater than the backward ones in QED. in QCD and in EW theory. 

Besides these 2 — 2, i.e. elastic processes, it is interesting also to study the 2 — > 2+71 -exclusive processes accounting 
Oh! for emission of n bosons accompanying the elastic 2 — + 2 annihilation. The point is that besides the conventional, 
^ ■ (soft) bremsstrahlung there can be emitted harder bosons. Emission of such bosons can be also studied in DLA to all 
orders in the couplings, providing the hard bosons are emitted in cones with opening angles <^ 1 around the initial 
e~^e~ beams, i.e. in the multi-Regge kinematics. In this case, the most important part of the inelastic scattering 
amplitudes accounting for emission of n bosons consists of the kinematic factor ~ (1/fci _L)---(l/fcn j_) multiplied by 
some function M which is called the multi-Regge amplitude of the process. The energy dependence of M is controlled 
I hy n + 1 electroweak Reggeons propagating in the crossing channel. Description of the multi-Regge photon production 
in the backward e+e~ — > ^i^ -annihilation was considered in Ref. |^ and in Ref. The multi-Regge amplitudes 
for giuon production in the backward annihilation of quark-airtiquark pairs were considered in Ref. | ^ . 

In the present paper we calculate the scattering amplitudes for electroweak boson production in e^e 
into quarks and leptons assuming that the bosons are emitted in the multi-Regge kinematics. We use the approach 
of Refs. and account for electroweak double- logaritmic contributions to all orders in the electroweak couplings. 

The paper is organised as follows: in Sect. II we consider emission of one EW boson in e^e~ -annihilation into a 
quark-antiquark pair. We compose the infrared evolution equations (IREE) for the amplitudes of these processes. 



annihilation 



The IREE are solved in Sect. Ill . A generalisation of these results to the case of emission of n bosons is given in 
Sect. IV . Emission of the EW bosons in e+e annihilation into leptons is considered in Sect. ^ . Results of numerical 



calculations are presented and discussed in Sect. ^ . Finally, Sect. VII is for conclusive remarks 
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II. EMISSION OF ONE ELECTROWEAK BOSON IN THE MULTI-REGGE KINEMATICS 



Let us start by considering the process e~^{p2)e^{pi) q{Pi) q{p'2) accompanied by emission of one electroweak 
boson with momentum k. Energies of the bosons are assumed to be ^ Mz- There are two kinematics for this process 
that yield DL radiative corrections. First of them is the kinematics where p'l ~ Pi, p'2 ~ P2- Obviously, 



[pi +P2f > h 



tl=q!^ ip'l-pif 



{P2 



p'2? 



(1) 



in this region. Eq. (|^) means that the final particles are in cones with opening angles 6 around the e^e beams. 
The second kinematics is the one where p'^ ~ p2 , P2 ~ Pi and therefore 



/2 / / 

ui=q 1 = [P2 



Pi) 



U2 



q'l = {p2-p',f 



(2) 



9-^1 around the e+e 
the u -kinematics. Both 



Eq. (H) means that the final particles are also in cones with the cmf opening angles tt 
beams. Through this paper we call kinematics (|l|) the t -kinematics and the kinematics 
of them are of the Regge type and studying them is similar in many respects. 

Instead of directly calculating inelastic amplitudes A'^^''^''^^ describing emission of any of 7, Z, W , it is possible 
to calculate first the amplitudes A*^"^ and A'^^^ (r = 1,2,3) describing emission of the isoscalar and the isovector 
bosons respectively. When expressions for such amplitudes are obtained, the standard relations between the fields 
7, Z, W and the fields corresponding to the unbroken SU{2) (g) U (1) can be used in order to express in terms 

of Aq, Ar- This way of calculating A^^'^''^^ is technically simpler than the direct one because when the radiative 
corrections are taken into account in DLA, contributions proportional to masses in propagators of all virtual EW 
bosons are neglected and therefore both the isoscalar and the isovector fields act as independent ones. It makes more 
convenient operating with virtual isoscalar and isovector bosons than with 7, 2', W -bosons. 

It is also convenient to discuss a more general process where lepton l^{pi) (instead of e~) and its antiparticle li'{p2) 
(instead of e"*") annihilate into the quark-antiquark pair q^{p[) qj'ip'2) and a boson. The emitted boson can be either 
the isoscalar boson A^ with c = or an isovector one A^ with c = 1, 2, 3. We consider first the most difficult case when 
both and are left-handed particles, transitions to the other chiralities are easy to do. The scattering amplitude 
of this process is qjq^ {M'^Y-Mi'l^ where the matrix amphtude {M'^)\i is the object to calculate. In order to simphfy 



operating with the isospin matrix structure of {M'^)\^,, it is convenient to regard the process in the crossing channel, 
i.e. in the t -channel for kinematics (|l|) and in the u -channel for kinematics (||). When the process IHii q^qjiAc is 
considered in the t -channel, its amplitude can be expressed through the same matrix (M^)*^^ , however with initial 
(final) t-state being qjP{Acq^ k') : 



q^ k,{M%}q,l 



(3) 



We have extracted the kinematic factor 2/k\ in order to simplify the matching condition (|3^) we will use. The 
initial cross-channel state qjl^ in Eq. (|^) can be expanded into the sum of the isoscalar and of the isovector irreducible 
SU{2) representations: 



qjl 



qbl" 



(4) 



The same is true for the final q^ -pair. Therefore, (M"^) - ^ can be represented as the sum 



k=0 



(5) 



of the invariant amplitudes Mk (k = 0, 1, 2, 3, 4), each multiplied the projection operator {P^)\y corresponding to an 
irreducible SU{2) -representation, fc = 0, 1 correspond to emission of the isoscalar field and k = 2,3,4 correspond 
to the isovector fields emission. Then, the projection operator (-Pp)!/ describes the case (see Fig. 1) when both the 
initial t -channel fermion state and the final one are SU (2) singlets. 

Obviously, in this case the emitted boson can be isoscalar only, i.e. c = 0. Therefore 



(6) 
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Mo Ml M2 



FIG. 1: The multi-Regge invariant amplitudes Mr (and the projector operators) in kinematics (1). The dotted lines correspond 
to the isoscalar Reggeons; whereas the zigzag lines stand for the isovector ones. The dashed lines denote isoscalar vector bosons 
and the waved line correspond to the isovector boson. 



The projection operators Pk, with {k = 1,2) describe the cases when both the initial and the final t -channel states 
are the isovector SU{2) states. However, Pi corresponds to the case when the emitted boson is isoscalar. 



iP^)lf,=2S^0itm)iit.n))- 

whereas P2 describes emission of isovector fields: 

(^):/-(4)^(n6a(ta)> 



(7) 



(8) 



T'^ (c = 1,2,3) in Eq. (||) stands for SU{2) generators in the adjoint (vector) representation. 

Projector P3 correspond to the case when the initial fermion state is the SU{2) singlet whereas the final one is 
the SU{2) vector. Projector P4 describes the opposite situation. The emitted boson is isovector in both these cases. 
Therefore, 



All operators in Eqs. (^, ^) are orthogonal: 



{PAYdiPBYd 



(9) 



(10) 



/ \2 



Below (see Eq. (^4|)) we will show that the invariant amplitudes M^jM^ do not have DL contributions. It leaves 
us with amplitudes Mo,i,2 to calculate. These invariant amplitudes account for radiative corrections to all powers in 
the EW couplings in the DLA. The arguments of Mk are 

si ^ {p[ + w 2pik , ti = 

S2 {p'2 + fc)^ ~ 2p2k , t2 

so that 

Si S2 = s/c^ 

The kinematics is the t -channel multi-Regge kinematics when 



1!- 



= (P2-P2)^ 



(11) 

(12) 
(13) 

Similarly, in order to simplify studying the isotopic structure of {M'^yd of Eq. (^ in kinematics (^, it is convenient 
to consider it in the u -channel where it can be expressed through u -channel invariant amplitudes M^: 



si 2 > <i,2 > Mi 



(14) 



k=0 



Operators {Pk')ij'i'j describe irreducible SU{2) representations, which for this channel are either symmetrical or 
antisymmetrical two-quark states: 



P ' 

-in 



51 5j, 



515], 



Pi 



5- 5i+5^5', 



5l{t% + (r)l5^,+5],(t^)i + (f)),5i 



(15) 
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Operators Pq, P{ describe emission of the isoscalar field whereas describes the isovector field emission. Let us 
note that operator 



^3 - 2 



st{^)], + (f)tsj,^s^,{ni-(f)pi (16) 



must be accounted for^ (cf Ref. 0). However, the invariant amplitude Afg related to P' does not yield DL contributions 
in the case of SU{2) though it does in the case of SU{N) with N > 2. It leaves us with three invariant amplitudes 
Mq I 2 (just like it was in the case of the t -kinematics). They depend on si.2 and on ui^2- In the multi-Regge 
kinematics (||) 

Sl,2 > Ui^2 > Ml ■ (17) 
In order to specify the multi-Regge t [u) -kinematics completely, we assume that 

> ^2 , {ui > U2). (18) 

The opposite case can be considered similarly. The kinematics where ti ^ t2 (wi ^ 1*2) means emission of soft 
electroweak bosons. This kinematics will be considered below separately. 

From the point of view of the Regge theory, accounting for radiative corrections in kinematics ([l^ , |l7| ) can be 
expressed through exchange of Reggeons propagating in the cross channels. Therefore operators Pqi (^o) of Eq. (||) 
(Eq. ([isl)) imply that amplitude Mq (Af/,) is controlled by two isoscalar Reggeons whereas the projection operators Pi 2 
of Eqs. (0,||) (operators P[ 2 of Eqs. (p^)) imply that the energy dependence of amplitudes Afi_2(Afi 2) is controlled 
by two isovector Reggeons. In contrast to it, one of the Reggeons in amplitudes Af3.4 is isoscalar and the other is 
isovector. 

Besides si^2 and ti^2 ("1.2), invariant amplitudes Afo,i,2 (-^^0 1 2) depend also on the infrared (IR) cut-off /i introduced 
in order to avoid IR singularities from integrating over virtual particle momenta. We use the IR cut-off /i in the 
transverse space. However, definition of /i for radiative amplitudes differs from the definition for elastic amplitudes. 
In this paper we introduce the same way as it was done in Let us denote kj°'\^ to be the component of a virtual 
particle momenta k[ transverse to the plane formed by momenta a and 6, with a ^ h. Then, the IR cut-off /i obeys 

M < (19) 

for alH = 1, . . . when a,b = pi,p2,p[,P2, k. In the present paper we assume that /i « Mz- 

In order to calculate Mr we generalise to the EW theory the technique applied earlier to investigation of the similar 
inelastic processes in QED[|| and in QCD§. The essence of the method is factorizing DL contributions from the 
virtual particles with minimal fc;°|^ and differentiating with respect to In/i^. At ti,t2 ^ /i'^, such particles can only 

be virtual EW bosons. Factorizing their contributions leads to the IREE for amplitude {M'^Y-j',. This equation is 
depicted in Fig. 2. 



, fe k k 




+sym 



FIG. 2: IREE for Mz- Letters inside the blobs stand for infrared cutoffs. 



^ The gluon fields In Ref. |7| were described in terms of matrixes (fA'^)^ 
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Applying to it the projector operators of Eqs. - 19[ |T5|) leads to the following IREE for Mr (see Refs. 
for technical details): 

dMr dMr dMr dMr 1 r, , / , 2> , . N 1 , . 

■ [hr Hs/^?) + hr{yi + 2/2) + TOfcJ/i] Mr . (20) 



2 



dpi dp2 dyi dy2 Stt 

We have used in Eq. ( |20| ) the fact that, according to our assumption Eq. (|l8|), k\ « ii(ui) and introduced the 
logarithmic variables 

Pi,2 = ln(si^2/Ai') , yi,2 = ln(ii,2/M') • (21) 
The numerical factors 6r, hr and in Eq. (kQ) are: 



g'^(y-n^ , , , , . 9'\Y-Y'f 
bo = , oi = 02 = 2.9 + , 

3g^ , g'^YY' , g'^YY' 

ho = — + , hi^h2^~- + , (22) 

rriQ = 'Til = , ra2 = g^ ■ 

The IREE for the invariant amplitudes Mq ^ 2 can be obtained similarly. It has the same structure as Eq. (|2^) , 
though everywhere ^1,2 should be replaced by U1.2. It means in particular that 7/1^2 should be replaced by y'l 2 = 
1^(^1.2)/ p,^- This replacement and replacement of operators -Po,i,2 by Pq ^ 2 results in replacement of the factors br, hr 
by b'r, h'r ■■ 



g'\Y + Y'f ^, ^ , , g'\Y^Y') 



l\2 



b'o = : &i - ^2 = 2.9 



4 

Factors mj, in the IREE for M'^ coincide with factors rrir in Eq. (p^). Therefore after replacements yi_2 — ^ 2j 
and &o,i,2 &o,i 2; /*o,i,2 ^0,1,2 -Eq. (E3) for amplitudes Mr in kinematics ([l3| ) holds for amplitudes describing 
e"'"e~ — > 99 -annihilation in kinematics (pT|). 



III. SOLVING THE EVOLUTION EQUATIONS FOR 



In order to solve Eq. (20), it is convenient to operate with the Mellin amplitude Fr related to through the 
Mellin transform: 

Mr = r ^^e-^''^+-^''^Fr{oJuu;2,yi,y2) • (24) 

In the Lu -representation, multiplying by pi corresponds to —d/duji. Using this and Eqs. ( |l^ , p^ ), we can rewrite 
Eq. (|^) as 



^dFr 


dFk\ 




\duji 


1 - 

duj2 ' 





-hr- mr)y\ - hry2] Fr . (25) 

For further simplification, it is convenient to introduce variables 2:1^2 and zi.2 : 

a;i,2 = , 21,2 = -Kyi,2 (26) 

where = y^br/Sn^ . 

In terms of Xi, Zi, the differential operator in the left hand side of Eq. ( p5| ) acquires symmetrical and simple form. 
Thus, we arrive at 

dFr dFr dFr dFr , o . Ir. 

^ + ^ + ^ + ^- [ixi+X2) + {l+Pr)zi+lrZ2]Fr (27) 
where (3r = —{hr + mr)/br, 7r = —hr/br- 
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The general solution to Eq. can be written as 

Fr = ^r{xi - 22, X2 - ^2, ^1 - ^2) GXp 



a;? + 



^+(l + /3.)^+7.f 



where unknown function has to be specified, ft can be done in particular through matching 

Fr{xi,X2,Zi,Z2)\z2=0 = Fr{xi,X2,Zi) 



(28) 



(29) 



where Fr is related through the Mellin transform (^4|) to amplitude Mr of the same process in the kinematics 
Eqs. (|l|Jl|) though with M^- The IREE (|^ for Fr differs from the IREE of Eq. (|^ for Fr in the following 

two respects. First, there is no Z2 dependence in Eq. (|3^). Second, in contrast to Eq. (|27|), the IREE for Fr contains 
an additional term (that we denote dQr{x) / dx) (and will specify below in Eq. (^2[)) in the rhs: 



(30) 



dFr dFr dFr r, ^ . ^ dQr(x2)i ~ 



This new term corresponds to the situation when the particles with the minimal transverse momenta are the t2 
-channel virtual quark pair (see Fig. 3). 




Pi -P2 
FIG. 3: The soft fermion contribution to the IREE for Mz 



This contribution is /i -dependent only when t2 ~ ^J.^■ The intermediate two-particle state in Fig. 3 factorizes 
amplitude Mr into a convolution of the same amplitude and the elastic amplitude Er- The explicit expressions for 
the elastic electroweak amplitude Er were obtained in Ref. Q . The particular case where the produced particles were 
a right handed lepton and its antiparticle was studied in Ref. [D . For all cases, the Mellin amplitude fr is related to 
Er through the Mellin transform 

Er - / ^^^'fr{^) : (31) 



— lOO 



with p = ln(s/^ ), can be expressed through the Parabohc cyhnder functions Dp^\x) with different values of Pr'- 

1 dMe^"/^Dj,M)] Dp.-i{x) 1 dQr{x) 

fr[x) = — = i: , , = . (32) 

Pr dx Dp^\x) Pr dx 

The term dQr{x2)/dx2 in the rhs of Eq. ( |30| ) corresponds to the contribution of the right blob in Fig. 3 to the IREE 
The general solution to Eq. (t50|) is 



Fr = ^r{xi - Zi,X2 - Zi) —— ^ CXp 

Dp^(x2)exp(a;^/4) 



't'^ I 

1 I 2 



2 ■ ^ ' ^'^ 2 



(33) 



where there is, again, an unknown function $r- 

In order to specify we use the factorisation of bremsstrahlung bosons with small k± which takes place (see 
Refs. §,§,§,0) both in Abelian and in non- Abelian field theories. In the context of the problem under consider- 
ation it states that when zi = 0, the radiative amplitude Mr and the elastic amplitude Er are related: 

MrU.^O^ErGr. (34) 

Er in Eq. (^J) are the invariant amplitudes of the elastic annihilation process (see [Q); Gr = g'{Y ± Y')/2 for r = 
(for invariant amplitude Mo(Mq) the sign is "-)-" ("-")); = 5 for r = 1,2,3. 
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Eq. (34) means that when zi = 0, the two Reggeons in every amphtude Mr converge into one Reggeon that controls 
Er energy dependence. However, such convergence is possible in the DLA only when both Reggeons are either isoscalar 
or isovector. This rules amplitudes M3.4 out of consideration. Obviously, this property of the multi-Regge amplitudes 
holds for the more complicated cases when the number of involved Reggeons is more than two. This property was 
first obtained in Rcf. and was called "Reggeon diagonality" . 

The matching (^) can be rewritten in terms of Mellin amplitudes (j)r{xiTX2) = Fr{uii^uj2i zi)\zj^=o and fr{(^) '■ 



/too ^ 
-too 2«7r 



friuj) . 



(35) 



We have used in Eq. ( |35| ) that according to Eq. (|12|) p — pi + p when zi = 0. For the amplitudes with positive 
signatures that we discuss in the present paper, the transform inverse to Eq. (|2^ ) can be written as 



Mr{pi,2,yi.2) 



—UJ1P1-UJ2P2 



Fr{uJi,UJ2,yi.2) 



(36) 



Applying this transform to Eq. (Ra) at zi = leads to 



/loo J <.oo 
— / dpidp2e^"""^)''^+^"-"^)''V.('^) = 



Gr 



too 
too 



duj 



2^7^ (cj — uJi){iLJ — UJ2) 



(37) 



Choosing the integration contours in Eq. (^7|) so that SRcj < 3fia;i^2 allows us to integrate over uj by closing the 
contour to the right, which does not involve dealing with singularities of fr. When this integration is done, we arrive 
at 



, /r(wi) - /r(w2) 
(pr = 

UJ2 - 

After that it is easy to obtain the following expression for M^- 

^ diOi dL02 I S2 ] 



(38) 



Mr = G 

Dp^{x2 - Zi) 



fr{xi) - fr{x2) 



LO2 — UJl 



exp 



(1-2/3,) 2 (I-27,) 



If we choose Ikxi < dix2, Eq. 



Dp^(x2 - Z2) 
J39I ) takes simpler form: 

Mj- — GrRr 



(39) 



(40) 



where 



Rr — 



duji duj2 ( si ^ 
2m 2m \q\ _ 



s-2 



1 -Pp,-i(a:i - zi) 

\Jqlql j uj2~uji Dp^ixi-zi) 



Dpr{x2 - zi) 
Dp^{x2 - Z2) 



■ exp 



(1-2/3,) 2 (1 - 27,) 2 

-2i . Z9 



4 ' 4 

The amplitudes M^'^\ M^^\ M^^^^ of the electroweak boson production are easily expressed through i?,; 
M^'') = cos ^vfM'"' + sin 6IvfM'^^ 

— g cos (^Rq + Ri ) + g sin 9wR2 , 



(41) 



sm tlw 



gsm6w{Ro + Ri) + 9 cos6wR2 , 



M'^^^^ = (1/\/2)[m(i) ±iM(2)] = (5/V2)i?2 (42) 

when the boson are produced in kinematics ( p^ ) . For kinematics ( p^ ) , the boson production amplitudes are expressed 
through R'^: 

m(^) = g cos 9w {Rq + R'l) + 9 sin 0wR'2 : 
M(^) = -gsmew{Ro + R'i) + 9COsewR'2 , 
M(^*) = {g/V2)R', . (43) 
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The exponent in Eq. ( [4l| ) is the Sudakov form factor for this process. It accumulates the soft DL contributions, 
with virtuaUties < zf. The harder DL contributions are accounted through Dp -functions. It is convenient to perform 
integration over uJi,2 by taking residues. Such residues are actually zeros Xk{r) (fc = 1, . . . ) of involved Dp^ -functions, 
so 

A x^(r) / \ ^r-Xk{r) 



'-h^> [7m) ■ 



Position of Xk{r) depend on values of pr in such a way that the greater pr, the greater are ^Xk{r). In particular, 
the real part of the rightmost zero = x{r) is positive when pr > 1. In other words, Rr increase with the total energy 
when Pr > I. Ref. states that 



, _ 3-YY't&ii'^0w , _ 1 + yy tan^ Qw 

^° " (y + r')'tan2 0w' Pi - ^2 - -g-f (y + rO'tan2 0w ^ ' 



and 



3 + tan^ 6lvi/ l-rr'tan2 6'vF 

^"^(T^YoW^' ^^=^^ = "8 + (r-r')^tan2^H^- 

Therefore, only Afo and A/q grow with increase of the annihilation energy whereas the amplitudes M\^2 and M' 
are falling. 

Let us discuss the asymptotics of i?r first. The asymptotics of the energy dependence of each is controlled by 
two identical isoscalar (isovector) leading Reggeons. Intercepts Aj (j = 5*, V, 5", V) of these Reggeons are related to 
the position of the rightmost zero x{j) of the Dp^ -functions so that 

As = AoS(po) , As' - y^xij,'^) , Ay = Aix(pi) , Ay, = A'ii(p'i) . (47) 



We remind that \j = ^Jh^JWi? . Therefore we arrive at the following asymptotics: 

Mo ~ g's^^ , Ml ~ g's^^ , ~ gs^^ , 

M^ ^ g's^s' , M[ ~ , ^ gs^v . (48) 

As the intercepts of the isoscalar Reggeons are greater than the ones of the isovector Reggeons, the asymptotics of 
the exclusive cross sections u'''^ and a'^^^ of the photon and Z -production is given by contributions of the isoscalar 
Reggeons with intercepts A5 and A5/. Therefore, the only difference between these cross sections is the different 
couplings of these fields to the isoscalar Reggeons. So, we conclude (see Eqs. (|4^) , (^) ) that asymptotically 

^^^^ 9 , ^ 

tan^ Ow ■ (49) 



r(7) 



Accounting for contributions of other zeros, x*-''-' changes the values of cr^'^^ and ct^^-' but does not change Eq. (|4S|). 
In contrast to Eq. (^9|), the asymptotics of the ratio a^'^'^/a^^'' depends on s. The point is that the exclusive cross 
section tr^^^ of W -production involves the isovector Reggeons with smaller intercepts. So, asymptotically this cross 
section obeys the following relation: 

(7) 



However, contributions of other zeros of Dp -functions change this asymptotic relation. Results of numerical 



calculation of and a^'^\ and accounting for non-leading DLA contributions are discussed in Sect. 



IV. EMISSION OF n VECTOR BOSONS IN THE MULTI-REGGE KINEMATICS 

The arguments of the previous Sects, can be extended in the straightforward way to the case when the e^e^ 
-annihilation into quarks is accompanied by emission of n isoscalar or isovector bosons with momenta fci,...,fc„ in 
the multi-Regge kinematics. It is not difficult to generalise expressions of Eqs. (^[]|, |l5|) for projection operators to 
the case of the n boson emission and obtain new projector operators. First of all, let us note that all non-diagonal 
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projectors should be ruled out of consideration by the same reason as it was done in Sect. 1. Therefore, the invariant 
amplitudes of emission of n bosons involve n + 1 identical intermediate Reggeons. The isotopic quantum numbers of 
the Reggeons depend on the initial fermion state and on the isospin of the emitted bosons. If the initial fermion state 
is isoscalar (or antisymmetric), the same is true for all intermediate Reggeons and therefore only isoscalar bosons can 
be emitted in these cases. The projector operators for this case are again Pq i (Pq i) with trivial adding factors 6^0 
for every isoscalar boson. If the initial fermion state is isovector (or symmetrical), the emitted bosons can be both 
isoscalar or isovector gauge fields. Accounting for emission of the isoscalar bosons does not require any changes of the 
projectors. When r (r < n) isovector bosons ci, C2, .., Cr are emitted in the t -kinematics, the operators P2 of Eq. (H) 
should be replaced by 

(i^'"'"---''"^):/ =(^;>)Mr^'-)'>a.---(T^^^)a.a,(T^^)a.a,(tal)} ■ (51) 

A similar generalisation of operator P2 is also easy to obtain. The new invariant amplitudes Mj{M'-) corresponding 
to these operators depend on n + 1 variables Si : 

si = 2piki , S2 = 2fcifc2 , . . . , s„+i = 2knP2 (52) 

and on ti{ui) in the case of the t (u) - kinematics: 

U = <i , Ut = 9'? (53) 

where 

qi = p'l-pi , 92 = gi - fcl , ■ • ■ , Qn+l = Qn - kn = P2 ^ P2 

Q'i = P'2-Pl ' 92 = ^'l - fcl , ■ • ■ , q'n+l =q'n- kn = P2 ~ p\ ■ (54) 

The kinematics is the multi-Regge t -kinematics if 

s, > > Ml (55) 

and it is the multi-Regge u -kinematics if 

S^ > Uj > (56) 

with i, j = 1, . . . ,n + 1. In order to define these kinematics completely, one should fix relations between different ti 
(different Ui). In this paper we consider the simplest case of the monotonic ordering. We assume that 

ii > t2 > • ■ • > in+i (57) 

for the case of the multi-Regge t -kinematics and the similar monotonic ordering 

Ul > U2 > ■ • ■ > (58) 

for the case of the multi-Regge u -kinematics ^. Eqs. (|5^ , |5^ ) read that in the both kinematics 

si . . . s„+i ^ sk\ ^...kl^ ^ . (59) 

It is also convenient to introduce variables pi = ln(si//^^) and yi where yi = ln(ii//^^) for the forward kinematics 
{y\ — ln(iii//i^) for the case of the backward one). In these terms. The IREE for looks quite similar to Eq. ^ 



dpi "' dpn+i dyi " ■ dy 



n+l 



1 

'8^ 



bj ln(s//i^) - h-iiyi + y^+i) + ^miyi 

I 



Mf , (60) 



^ Scattering amplitudes for other multi-Regge kinematics can be calculated similarly (see Ref. [0). It is worth to mention that amplitudes 
for kinematics (pTl) and (Bq) yield main contributions to the inclusive cross section when integration over the EW boson momenta is 
performed. 
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where bj, hj are given by Eqs. ( |22|J23| ); mi = g if the boson I with momentum ki is isovector, otherwise m = 0. Let 
us consider for simphcity the case of emission of isoscalar bosons. Introducing the Melhn ampUtude through the 
transform 



Mn = 



2Tn 



F„(cJi, . . .,UJn+l,yi, ■ ■ ■,yn+i) 



and using notations Xi , Zi defined as 



(61) 



(62) 



we transform Eq. (61) into the following one: 

dFr, , , dF„ , dFn 



OF,, 



with the general solution 



dxi dxn+i dzi dzn+i 

[{Xi H h Xn+l) + Zi-\ \- Zn + h(zi + Z„+l)] F„ 



$n (a^l ~ Zn+1, X2 — ^n+l, • ■ ■ , Xn+1 — Zn+lJ ^1 — Zn+1, ■ ■ ■ , Zn — Zn+l) 



(63) 



exp 



Sn+l{x) + Sn{z) + h{- 



'1 ^ N 



where we have denoted Sr{a) = '^f /2 • An unknown function $„ can be specified through the matching 



(64) 



(65) 



where the Mellin amplitude F„ describes the same process in the multi-Regge kinematics (|5^,^,^) though with 



ql^^ = The IREE for F„ is 



OK 
dxi 



+ 



dF^ dFn 



dxn+i dzi 



dFn 

dz„+ 



(Xl H 1- Xn+l) + Zl-\ 1- Zn+1 + hzi + 



dQr{Xn+l) 



Fi 



(66) 



dXn+i 

where Q is defined by Eq. ( ^^ . The general solution to Eq. ( |6^ ) can be obtained quite similar to the one of Eq. (|6^) : 

Fji — (-^1 ^n: ^2 ^n; ■ ■ ■ i -^n+l -^■ni -^1 -^n? ■ ■ ■ 5 -^n — 1 ^n) 

exp[5„+i(a;) + 5„(z) + ;i^] . (67) 

It also contains an unknown function $. In order to specify it we use factorisation of the photons with small k±: 

M„U„=o = M„-i(si, . . . ,Sn;ql, ■ ■ ■ , ql_i,fi'^) ■ (68) 

Rewriting this equation in terms of the Mellin amplitudes and performing the transform inverse to Eq. (|6l]), we 
express Fn through amplitude Fn^i: 



1 



Fn—1 (^li ■ ■ ■ ;^n — 1; ^ri; -21 , . . . , Zfl—l) 



F„ 



Combining this equation with Eqs. ( |61|J67| ) leads to the following recurrent formula for Fn'. 

Fyi (xi , . . . , Xn-\-l j ^1 5 ■ ■ ■ 1 ) ~ Fn~ 1(^1 ; ■ ■ ■ ; i ^1 ; ■ ■ ■ ; ^n — l ^n) 



(69) 



Dp{Xn+l — 



exp 



Sn+lix) - Sn+lix - Zn) + Sn{z) - Sn-liz - Zn) + ^[zf + - {zi - Z„)^] 



(70) 
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Using this formula leads to the following expression for the amplitude M^-"'' of emission of n isoscalar bosons in the 
ordered kinematics (pTUs 



Mi") = 



,{Y±Y') 



duji 



dLJn+: 



exp 



Sn+l 



[hj. — 2hr) 
46l 



(zf + zj) 



(71) 



Dp^{x2 — Z2) Dp^{Xn^i — Z„+i) 

where Y zLY' corresponds to the kinematics (^,^8|) respectively. Eq. ( [7l|) implies that the contours of integrations 
obey SRxi < .. < 3fix„+i. After that one can perform integration in Eq. ( |7l| ) by taking residues in the Dp^ zeros. 

When k of the isoscalar bosons are replaced by the isovector ones, {g'{Y ± Y')/2)" in Eq. ( ^l| ) should be replaced 
by g''{g'{Y zL Y')/2)^~^; the factor {in/2br)zf for each of the emitted isovector bosons should be added to the last 
exponent. Using the standard relation between gauge fields Ar and 7, Z, W, one can easily rewrite the gauge boson 
production amplitudes of Eq. ( [zil) in terms of amplitudes for the electroweak bosons production. Asymptotics of the 
scattering amplitudes of the photon and Z -production are governed by the isoscalar Reggeons whereas W -production 
involves the isovector Reggeons. Eq. (|7^) can be used for obtaining different relation between cross sections of different 
radiative processes in the multi-Regge kinematics. For example, 



(nZ) 



tan^" 



w 



whereas the energy dependence of the ratio cr^"'*'' /cr^"*^) is less trivial. Asymptotically, 



-0.36 



JnW) 



(72) 



(73) 



Results of accounting for the non-leading Reggeon contributions for (j("''''/cr'^"'^) can be obtained from Fig. 6 



because cr(«7)/cr(«M^) ^ {eV2/gy'-^a^'^'> /a^^l In obtaining Eqs. (@|7|) from Eq. ^ we have used that according 



to Eq. (M), (si) 



EMISSION OF EW BOSONS IN e+e -ANNIHILATION INTO LEPTONS 



Inelastic annihilation of e+e , with the e being left handed, into another left handed lepton I' (i.e. into fi or t) 
and its antiparticle I' can be considered quite similarly to the annihilation into quarks discussed above. In particular, 
explicit expressions for the new invariant amplitudes, Lr"'' of this process can be obtained from Eqs. (|39|) , (pO[) , ([7l|) 



by putting Y' = Y = —1 in Eqs. (p^),(p3[) for the factors br, 6^, h'^. However having done it, we obtain that 60 = 
(see Eq. (p^). It means that the IR evolution equations for the scattering amplitude Lq^^ of the inelastic annihilation 
e"'"e~ — > I'l' + ni7 + {n — ni)Z in the kinematics ([^) do not contain contributions proportional to ln(s///^) in the rhs 
and therefore the Mellin amplitudes /q""* (related to Lq through the Mellin transform (p4[), Eq. (|74|)) do not have the 
partial derivatives with respect to ujj (cf Eq. ([7l|)). In order to obtain expressions for the new scattering amplitudes 
Lq"\ let us consider first the simple case of emission of one isoscalar boson accompanying the forward e+e~ - 
-annihilation, assuming that both e~ and I' are left particles. It is obvious that for this case, the IREE of Eq. 
for scattering amplitude Mq has to be replaced by the simpler one, 



dL 



(1) 



dL, 



(1) 



dL 



(1) 



dL 



(1) 



dpi 



where we have denoted /iq — (3^^ 
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dp2 
")/4 . 



dyi 



dy2 



1 

'8^ 



ho{yi + y2)io 



(1) 



In terms of the Mellin amplitude /i^'*, Eq 



(wi 



|) takes the following form: 



dyi dy2 



= hoiyi + y2)fo 



I'V 

(I20I) 



(74) 



(75) 



^ the kinematic factor is also extracted from A^"' like it was done for amplitudes A4^"K 
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The solution to Eq. (|7q) respecting the matching condition (EJ) is (cf Eq 



(1) _ 9'Y f"^ dui du2 (siYUS2Y^ - h{x2)) 



2i7r 2i7r Vgj 



mm 



(tJ2 - t^l) 



„(l/87r2)/o(t^2)(yi-!y2)„-('»o/2)(y?+y^) 



where the Mehin amphtude /o = for the elastic e'^e —>■ i-i^ annihilation is (see 

/o = 47r2 L - X2-(3.g2+5'2)/8^2 



(76) 



(77) 



The last exponent in Eq. ( [76|) is the Sudakov form factor accumulating the DL contribution of the soft virtual EW 
bosons only. The other terms in the integrand account for harder contributions. The leading singularity (intercept), 
ujQ of the integrand of Eq. (fz^) is given by the position of the branch point of the rhs of Eq. (77). Therefore we obtain 



Wo 



27r Vsin^e 



1 



w 



cos^ 9w 



= 0.13 , 



so asymptotically 



„0.13 



(78) 



(79) 



The invariant amplitudes Lq"'' for production of n isoscalar bosons in the kinematics ( p7| ) when e^e annihilate 

(n) 

into another lepton pair can be obtained similarly. The IREE for the amplitudes Lq is 



dL, 



(«) 



in) 



dpi 

and its solution is 

4"'^(^)" 



dp, 



•n+l 



in) 



dyi 



+ .. 



dL 



in) 



dy 



71+ 1 



2^7^■ 



- yn+l)LQ 



/o(wi) 



in) 



(80) 



(W„+1 - W„)...(W2 - t^l) 



2i7r \q(. 

=,(i/8^^)[/o('^2)(ai-a2)+...+/o('^„+i)(s/„-y„+i)]g-(So/2)(j/?+..+af.+i) 



n) 



when 5ft(wi) < 5R(wi+i), « = 1, .., n . Their asymptotic s-dependence is also given by Eq. ([79[). The results of numerical 
calculations for the cross section of 7, Z and W production in e+e^ I'V are presented in Fig. 6. 



VI. NUMERICAL RESULTS 



In order to estimate at what energy scale one might hope to observe the predicted asymptotical behaviour of cross 
sections of exclusive and Z, 7 production we have first to account for all non-leading DLA amplitudes for left 
and right chiralities of initial e^e~ and final qq or II pairs. There are many such amplitudes, but all of them can be 
easily calculated as described in previous sections. The results for Regge intercepts for the forward (i-channel) and 
backward (u-channel) kinematics are collected in Table | for the final qq and in Table || for the final II. 

Evidently, in far asymptotics the leading contribution for production comes from Fi (isotriplet) of the backward 
eje£ qhqL whereas the leading contribution to {Z, 7) production comes from i^o (isosinglet) of the forward eje£ — > 
qLqL- However differences between the non- leading and leading intercepts are small, and one can expect the role of 
the first to be essential at real energies scales. Moreover, the effects of non-leading intercepts of the same amplitude 
can be also large enough at real energies. Therefore it seems reasonable to numerically compute the energy dependent 
amplitudes, Mr, by taking the inverse transform of the IREE solutions Fr{uj), and to calculate with them the inelastic 
cross sections of boson production in central region (in cmf) with p? — M|. It seems also suitable to sum over 

the final qq or II isotopic states, fixing only the emitted boson isotopic state. 

Easy but cumbersome calculations lead to the following results. production in the forward and backward 

kinematics is described by the same formula (though with different amplitudes involved!): 

^-^ = cf[\MLLT? + \MnL?] , (82) 
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TABLE I: Rightmost zeros xq of parabolic cylinder functions Dp{x) determining the values of the leading singularities uiq of 
different Mellin transform amplitudes Fr{uj) for e^e~ qq annihilation in forward and backward kinematics. 





forward kinematics 


backward kinematics 


V 




A 




P 




A 




-t'LLS 


5.796 


3.23 


0.026 


0.083 


24.68 


8.65 


0.013 


0.111 


Fllt 


-0.129 


-2.52±1.62i 


0.106 


-0.267±0.171i 


-0.805 


-1.98±2.62j 


0.039 


-0.076±0.101j 


FlRu 


-0.083 


-2.65±1.48i 


0.077 


-0.205±0.115i 


0.124 


-1.85 


0.063 


-0.117 


FhRd 


0.062 


-2.25 


0.063 


-0.142 


-0.05 


-2.81±1.35i 


0.071 


-0.199±0.095i 


Frl 


-0.042 


-2.87±1.31i 


0.077 


-0.222±0.102j 


0.05 


-2.365 


0.071 


-0.167 


FrRu 


-0.24 


-2.33±1.86i 


0.064 


-0.15 ±0.12 I 


6. 


3.32 


0.013 


0.0428 


FRRd 


0.75 


-0.34 


0.026 


-0.009 


-0.188 


-2.40±1.76i 


0.051 


-0.124±0.090j 



TABLE IL Rightmost zeros xq of parabolic cylinder functions Dp [x) determining the values of the leading singularities u)q of 
different Mellin transform amplitudes Fr{uj) for e'^e^ — > annihilation in forward and backward kinematics. Notations for 
isodublet components of I are taken as for muon doublet. 



Fr{uj) 


forward kinematics 


backward kinematics 


P 


Xo 


A 




P 


Xo 


A 




Flls 


oo 






0.132 


2.41 


1.32 


0.039 


0.051 


Fllt 


-0.090 


-2.63±1.50i 


0.103 


-0.270±0.154j 


-0.602 


-2.06±2.40j 


0.053 


-0.109±0.127j 


Flru 


















FlR^i 


0.172 


-1.64 


0.066 


-0.108 


-0.102 


-2.59±1.54i 


0.086 


-0.221±0.132i 


Frl 


0.172 


-1.64 


0.066 


-0.108 


-0.102 


-2.59±1.54 2 


0.086 


-0.221±0.132z 


FrRv 


















Frr^ 


oo 






0.077 


-0.25 


-2.32±1.88i 


0.077 


-0.179±0.145i 



where cto is the common Born cross section of the elastic process (see Mllt denotes M2 ampUtude of eje£ iLflL 
and MfiL denotes the amphtude of e^e^ — > (and similar for annihilation to leptons). Let us remind that 

are produced first as isovector boson Ai, A2 states and then transform to observable boson states. In contrast, {Z, 7) 
are being produced first as isoscalar B or isovector fields, and then transform to the observable states, Z mainly 
comes from A3 and 7 - from B. Easy but cumbersome calculation lead to the following cross sections for production 
of B and A3 bosons: 

a{Z,j)^a{A3) + a{B) , (83) 



(84) 



where again the amplitudes M involved are either forward or backward amplitudes: Mllt = M2 of eje^ QlQl, 
Mlru stands for eje^ urur, MLRd stands for eje£ dRdR and Mrl for e^e^ qLQL, and 



1, 



\ I MLR.U I 



(■^Ci? "F^uj?)^ I , , |2 I (^en T^dn)^ I |2 

-. VV1rru\ H \^v^RRd\ 



(85) 



where "=f" denotes that "-" sign corresponds to forward amplitudes and sign to backward amplitudes denoted 
above. 

The same formulae can be used for e~ II annihilation channel: one has to substitute the appropriate amplitudes 
Mr and to replace electro- weak charge Yq with the appropriate Yi . 

As the Regge kinematics is dominating in the cross sections, we sum the contributions of forward and backward 
kinematics in what follows. The results of numerical calculations presented in Figs. ^, || show that at energies 
^/s < 10^ GeV and {Z,'y) are mainly produced in e+e" annihilation. And only at ^/s > 10^ 10^ GeV 

their yields from e'^e~ qq annihilation become greater (see Fig. m. 
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10^ 10" 10^ 10^ 10 
^/S /GeV/ 



FIG. 4: Dependence of exclusive and (^,7) production on the total energy of e'''e annihilation. The cross sections are 
divided by the differential elastic Born cross section ctq to make differences in energy dependencies more clear. 



(T 




10^ 10* 10^ 10*^ 10 



75 /GeV/ 



FIG. 5: Total energy dependence of W and (2', 7) production in different channels of e'^e annihilation: e'^e II - solid 
curves and e'^e~ ^ qq - dashed curves. 



The explicit asymptotical dominance of exclusive channel qq + (Z, 7) over the channel II + {Z, 7) stems from the 
fact that despite the leading Flls amplitude in the Table || has the greater intercept loq « 0.132 than ujo w 0.111 of 
the leading Flls in the Table ||, its contribution is multiplied by the zero factor in Eq. (|85|). 

The numerical calculation for the ratio of to {Z, 7) production summed over the both annihilation channels 
e+e^ n and e+e^ qq is shovifn in Fig. ^. Thus the DLA predicts rather slow energy dependence of the ratio till 
y/s ~ lO^GeV and then its relatively rapid decrease. 

The energy dependence of the ratio Z to 7, 

is shown in Fig. |^. In far asymptotics radiation of isoscalar field B dominates over radiation of isovector field A and 
the ratio tends to the fixed value tan^ 9w ~ 0.28. 
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'^(2,7) 
. 8r 

0.7- 
0.6- _ 
0.5- 
0.4- 
0.3- 
0.2- 
0.1- 

10^ 10* 10^ lO*' 10 
/GeV/ 



FIG. 6: Total energy dependence of to (^,7) rate in e^e annihilation. 

1 

. 8 
. 6 
. 4 



0.2 



10^ 10" 10^ 10^ 10^ 
VS /GeV/ 



FIG. 7: Total energy dependence of Z to y rate in e'^e annihilation. The dashed line shows the asymptotical value of the 
ratio: tan^ 6w ~ 0.28 . 

Worthwhile to note that apart from the results of j^] for pure QED, the figures of the Table |l^ show that backward 
eje£ — > IlIl amphtude (i.e. when antilepton foUows the direction of initial electron) for isoscalar channel in EW 
theory has the positive intercept though small enough if compared to forward annihilation amplitude. 

Let us emphasise that the demonstrated in Fig. excess of Z production over 7 production at y/s < 10'^ 10^ GeV 
as well as excess of {Z, 7) emission over emission in the same energy range shown in Fig. ^ and the dominance 
of II channel over qq channel shown in Fig. || may all happen to be just artifacts of the DLA. To get more reliable 
predictions for the cross sections one has at least to account for single logarithmic corrections as well. The presented 
figures show that account of the non-leading DLA effects can make observation of the theoretically correct predictions 
( p9| ) and (|5^) hardly possible even in far future. 
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VII. CONCLUSION 



In the present paper we have obtained exphcit expressions for the scattering amphtudes for the e+e^ annihilation 
into quarks and into leptons at the annihilation energies -^/i ^ 100 Gev accompanied by emission of n electroweak 
bosons in the multi-Regge kinematics, i.e. in the kinematics where the final particles are in cones with opening angles 
^ 1 around the initial e+e" beams. We accounted for the double-logarithmic contributions to this process to all 
orders in the EW couplings. We have shown that it is convenient to calculate amplitudes of this process in terms 
of the isoscalar and of the isovector amplitudes. The isoscalar amplitudes describe production of the isoscalar gauge 
fields. They are controlled by n + 1 isoscalar Reggeons propagating in the crossing channel. The leading intercepts 
of these Reggeons are positive (A5/ = 0.11 and A5 = 0.08) and therefore such scattering amplitudes grow when s 
increases. The isovector amplitudes bring sub-leading contributions to the production of the isoscalar bosons and 
in the same time give the leading contributions to production of the isovector gauge fields. They are governed by 
n + 1 isovector Reggeons with negative intercepts Ay/ = —0.08 and Ay = —0.27. It means that the amplitudes for 
isovector production decreases when s grows. These results lead in particular to the fact that production of each Z 
boson is always accompanied by production of a hard photon with the same energy ^ 100 Gev. In DLA, such hard 
photons are never produced without Z bosons. The cross sections of production of these photons and the Z bosons 
have identical energy dependence, however they are different numerically due to difference in the couplings. They are 
related by Eq. ( [72| ) at asymptotically high energies (> lO"^ Gev). The s -dependence of the ratio (T"^/tT"^ for lower 
energies is given in Fig. 0. The energy dependence of cross section for the W production is weaker than the one for the 
photons and the Z bosons by factor g^f^ asymptotically high energies. The s -dependence of these cross sections 

is shown in Figs. Through this paper we consider only the monotonically ordered multi-Regge kinematics js^ ) 
and (|5^). Accounting for the other kinematics can be done in a similar way. Though it is likely to bring corrections 
to explicit formulae for the invariant amplitudes Mr"\ it cannot change the asymptotic relation of Eq. (72) and the 
fact (see Eq. (f^) that cr(«7)/o-(«w') decreases with s. 
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